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Abstract  -  We  propose  a  large-scale  quantum  computer  architecture  by  more  easily  stabilizing 
a  single  large  linear  ion  chain  in  a  very  simple  trap  geometry.  By  confining  ions  in  an  anharmonic 
linear  trap  with  nearly  uniform  spacing  between  ions,  we  show  that  high-fidelity  quantum  gates 
can  be  realized  in  large  linear  ion  crystals  under  the  Doppler  temperature  based  on  coupling  to  a 
near-continuum  of  transverse  motional  modes  with  simple  shaped  laser  pulses. 
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Trapped  atomic  ions  remain  one  of  the  most  attractive 
candidates  for  the  realization  of  a  quantum  computer, 
owing  to  their  long-lived  internal  qubit  coherence  and 
strong  laser-mediated  Coulomb  interaction  [1-4].  Various 
quantum  gate  protocols  have  been  proposed  [1,5-9]  and 
many  have  been  demonstrated  with  small  numbers  of 
ions  [4,10-14].  The  central  challenge  now  is  to  scale  up  the 
number  of  trapped  ion  qubits  to  a  level  where  the  quantum 
behavior  of  the  system  cannot  be  efficiently  modeled 
through  classical  means  [4].  The  linear  rf  (Paul)  trap  has 
been  the  workhorse  for  ion  trap  quantum  computing,  with 
atomic  ions  laser-cooled  and  confined  in  ID  crystals  [1-4] 
(although  there  are  proposals  for  the  use  of  2D  crystals  in 
a  Penning  trap  [15]  or  array  of  microtraps  [6]).  However, 
scaling  the  linear  ion  trap  to  interesting  numbers  of  ions 
poses  significant  difficulties  [2,4].  As  more  ions  are  added 
to  a  harmonic  axial  potential,  a  structural  instability 
causes  the  linear  chain  to  buckle  near  the  middle  into 
a  zigzag  shape  [16,17],  and  the  resulting  low-frequency 
transverse  modes  and  the  off-axis  rf  micromotion  of  the 
ions  makes  gate  operation  unreliable  and  noisy.  Even  in 
a  linear  chain,  the  complex  motional  mode  spectrum  of 
many  ions  makes  it  difficult  to  resolve  individual  modes 
for  quantum  gate  operations,  and  to  sufficiently  laser  cool 
many  low-frequency  modes.  One  promising  approach  is 
to  operate  with  small  linear  ion  chains  and  multiplex 
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the  system  by  shuttling  ions  between  multiple  chains 
through  a  maze  of  trapping  zones,  but  this  requires 
complicated  electrode  structures  and  exquisite  control  of 
ion  trajectories  [2,18]. 

In  this  paper,  we  propose  a  new  approach  to  ion 
quantum  computation  in  a  large  linear  architecture  that 
circumvents  the  above  difficulties.  This  scheme  is  based 
on  several  ideas.  First,  an  anharmonic  axial  trap  provided 
by  static  electrode  potentials  can  stabilize  a  single  linear 
crystal  containing  a  large  number  of  ions.  Second,  tightly 
confined  and  closely  spaced  transverse  phonon  modes 
can  mediate  quantum  gate  operations  in  a  large  archi¬ 
tecture  [19],  while  eliminating  the  need  for  single-mode 
resolution  and  multimode  sideband  cooling.  Third,  gate 
operations  on  the  large  ion  array  exploit  the  local  char¬ 
acter  of  the  laser-induced  dipole  interaction  that  is  domi¬ 
nated  by  nearby  ions  only.  As  a  result,  the  complexity  of 
the  quantum  gate  does  not  increase  with  the  size  of  the 
system,  and  the  gates  can  be  performed  in  parallel  on  ions 
in  different  locations  of  the  chain. 

The  proposed  ion  architecture  is  illustrated  in  fig.  1. 
It  is  a  large  linear  array  where  the  strong  confinement 
in  the  transverse  ( x,y )  direction  is  provided  by  the 
ponderomotive  Paul  trap  with  an  effective  potential  of  the 
form  V(x,y)  =  (m(jj2x2  +  mu)2y2) /2,  where  m  is  the  mass 
of  each  ion.  The  ions  are  initially  Doppler  cooled,  with 
a  number  of  ions  at  the  edges  of  the  chain  continuously 
Doppler  cooled  in  order  to  overwhelm  any  heating  that 
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Fig.  1:  (Color  online)  Linear  architecture  for  large-scale  quan¬ 
tum  computation,  where  lasers  address  individual  ions  and 
couple  to  local  modes  of  the  ions,  while  the  edge  ions  are 
continuously  Doppler  laser-cooled.  In  a  large  ion  chain,  more 
efficient  sympathetic  cooling  could  be  achieved  with  cooling 
ions  sparsely  distributed  in  the  ion  chain. 


occurs  during  the  gate  operation.  The  middle  portion 
and  majority  of  the  ion  chain  is  used  for  quantum 
computation.  Given  an  appropriate  axial  static  potential 
V(z)  from  the  axially  segmented  electrodes,  we  assume 
these  computational  ions  are  distributed  nearly  uniformly, 
with  a  neighboring  distance  of  about  ^  10  /xm.  This 
enables  efficient  spatial  addressing  with  focused  laser 
beams  along  the  transverse  direction  for  quantum  gate 
operations. 

When  the  axial  potential  takes  the  conventional 
harmonic  form  V(z)  mmuj2z2  / 2,  the  ion  array  is  subject 
to  the  well-known  zigzag  transition  unless  the  trap 
anisotropy  is  at  least  uox,y/uz  >  0.77N/y/logN  [2,17,20], 
where  N  is  the  number  of  ions.  As  N  becomes  large,  this 
structural  instability  occurs  first  at  the  minimal  distance 
at  the  trap  center  due  to  the  spatial  inhomogeneity  of 
the  ion  distribution.  When  the  ions  under  a  similar  trap 
anisotropy  are  instead  uniformly  spaced  by  neighboring 
distance  do,  it  is  easy  to  see  that  the  linear  structure 
is  always  stable  even  for  an  infinite  chain  so  long  as 
>  7C(3)e2/(2radjj)  «  4.2e2/(radjj),  where  ((l)  is  the 
Reimann  zeta  function  and  e  is  the  charge  of  each  ion. 
Therefore,  a  large  linear  structure  can  be  more  easily 
stabilized  so  long  as  static  potentials  from  the  trap 
electrodes  are  designed  to  accommodate  equally  spaced 
ions.  The  uniformity  of  the  distribution  is  critical  for 
scaling  because  the  minimum  ion  spacing  can  be  kept 
constant  to  avoid  the  zigzag  transition,  while  at  the 
same  time  the  maximum  ion  spacing  does  not  grow,  as  is 
required  for  operating  entangling  gates  efficiently  across 
the  whole  chain,  as  described  below. 

To  illustrate  the  general  method,  here  we  consider 
an  explicit  example  with  a  quartic  potential  V (z)  = 
a2Z2/2  +  a^z^/A  that  can  be  realized  with  a  simple  five- 
segment  electrode  geometry  as  shown  in  fig.  2(a).  Under  a 
quartic  trap  V(z),  the  axial  equilibrium  position  Zi  of  the 
i-th  ion  can  be  obtained  by  solving  the  force  balance 
equations  dU/dzi  =  0,  where  V  =  YUWizi)  +  y)\  + 

J2i<j  e2/\r i  —  Yj  |  is  the  overall  potential  including  the  ions’ 


Fig.  2:  (Color  online)  (a)  Sample  five-segment  linear  ion  trap 
with  voltages  Vi  (i  =  1,2)  to  produce  a  quartic  axial  potential. 
The  ions  are  confined  in  the  central  segment,  (b)  The  variance 
of  the  ion  spacings  sz  in  a  linear  quartic  trap  as  a  function  of 
the  trap  parameter  B  that  characterizes  the  ratio  of  quadratic 
to  positive-quartic  nature  of  the  potential,  (c)  The  distribution 
of  the  ion  spacing  at  the  optimum  value  B-—  6.1.  The 
computational  ions  are  within  the  dashed  lines,  where  the 
spacing  is  essentially  uniform. 

mutual  interactions.  We  optimize  the  dimensionless  ratio 
B  =  |o^2  / e2 1 2//3  (0^2  / <^4)  characterizing  the  axial  potential 
to  produce  a  nearly  uniformly  spaced  crystal.  To  be 
concrete,  we  consider  an  array  of  120  ions,  with  10  ions 
at  each  edge  continuously  laser  cooled  and  100  qubit 
ions  in  the  middle  for  coherent  quantum  gate  operation. 
We  solve  the  equilibrium  positions  of  all  the  ions  under 
V (z)  and  minimize  the  variance  in  ion  spacing  sz  = 

\/ Iho  Siiii  —  A zn)2  for  the  qubit  ions,  where  A zn 

is  the  distance  between  the  n-th  and  (n  +  l)th  ion  in  the 
chain  and  A zn  denotes  its  average.  The  variance  in  spacing 
is  shown  in  fig.  2(b)  as  a  function  of  the  parameter  B. 
The  value  of  sz  is  fairly  insensitive  to  B  and  reaches  a 
minimum  when  B  ^  —  6.1.  Here,  the  distribution  of  ion 
spacing  zn  is  shown  in  fig.  2(c),  which  is  remarkably 
homogeneous  for  the  qubit  ions  even  though  we  have 
optimized  just  one  control  parameter:  sz/Azn  deviates  by 
only  3%  over  the  entire  crystal.  In  this  configuration,  if 
we  take  Azn  =  10  fim  for  atomic  Yb+  ions,  we  only  need  a 
transverse  center-of-mass  frequency  (jOx^/2t:  >  221kHz  to 
stabilize  the  linear  structure.  In  this  paper,  we  actually 
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take  ujx/ 2tt  =  5  MHz,  as  is  typical  in  experiments  [21],  and 
such  transverse  confinement  would  be  able  to  stabilize 
linear  chains  with  thousands  of  ions  under  an  optimized 
quartic  potential.  We  note  that  the  ion  spacing  can  be 
made  even  more  uniform  by  adding  higher  order  multipole 
potentials. 

We  now  describe  quantum  gate  operations  with  this 
large  ion  chain,  mediated  by  many  transverse  phonon 
modes.  Given  the  equilibrium  positions  of  the  ions,  we 
can  efficiently  determine  all  axial  and  transverse  phonon 
modes.  We  then  apply  a  spin-dependent  laser  force,  with 
the  resulting  interaction  Hamiltonian  [22,23] 

#  =  FI  COS (kqn  +  fit),  (1) 

n 


where  the  transverse  displacement  qn  of  the  n-th  ion  in 
the  x-direction  is  expressed  in  terms  of  phonon  modes  a\~ 
with  eigenfrequency  Uk  and  the  normal  mode  matrix  b\ 
by  qn  =  Efc  bn^H/ 2 mwfe (a\ el“kt  +  ake~luJkt).  The  normal 
mode  matrix  6^  and  its  eigenfrequency  uok  are  determined 
by  solving  the  eigenequations  ^in^n  =  >  where 

Ain  =  d2U  / dxidxn  are  calculated  at  the  ions’  equilibrium 
positions  zi.  In  eq.  (1),  a *  is  the  Pauli  spin  operator  for 
the  n-th  ion,  Qn(t)  denotes  the  Rabi  frequency  of  the 
laser  pulse  on  the  n-th  ion  with  detuning  fi  from  the 
qubit  resonance,  and  the  effective  laser  momentum  kick 
k  is  assumed  to  be  along  the  transverse  x-direction.  (For 
twin-beam  stimulated  Raman  laser  forces  and  hyperfine 
state  qubits,  the  effective  laser  kick  carries  momentum 
along  the  difference  wave  vector  ki  —  k2  of  the  two  beams.) 
Due  to  the  strong  transverse  confinement,  the  Lamb-Dicke 
parameter  rjk  =  \k\^/h/2mu0k  1,  and  the  Hamiltonian 

H  can  be  expanded  as  H  =  -^jnk%xn{t)gt{a\e,'UJkt  + 
ake~iUkt)a zn  with  g\  =  gkb kn  and  x„(f)  =  Q,n(t)  sin  (fit)  (the 
effect  of  higher-order  terms  in  the  Lamb-Dicke  expansion 
will  be  estimated  later).  The  corresponding  evolution 
operator  is  given  by  [23] 


Fig.  3:  (Color  online)  (a)  The  gate  infidelity  SF  as  a  function  of 
the  laser  detuning  /jl  from  the  qubit  resonance,  with  optimized 
Rabi  frequencies  over  M  —  5  equal  segments  of  the  laser  pulse. 
Different  curves  correspond  to  different  gate  times.  The  120 
transverse  phonon  modes  are  all  distributed  within  the  narrow 
frequency  range  indicated  by  the  two  vertical  lines,  (b)  The 
shape  of  the  laser  pulse  r\Sl  (in  units  of  cjx)  that  achieves  the 
optimal  fidelity  (with  £F  =  8.5xlO-6),at  the  detuning  shown 
by  the  arrow  in  (a)  (with  (fi  —  ujx)/ujx  =  9.3  x  10-3)  and  the 
gate  time  t  =  500to.  The  dashed  (or  dotted)  lines  represent 
the  approximate  optimal  solutions  of  the  laser  shape  where 
only  4  (or  8)  ions  (from  the  59th  to  62nd  or  the  57th  to  64th, 
respectively)  are  allowed  to  vibrate  and  all  the  other  ions  are 
fixed  in  their  equilibrium  positions.  This  approximation  does 
not  significantly  change  the  optimal  laser  pulse  shape  compared 
to  that  of  the  exact  solution  (represented  by  the  solid  line) 
where  all  ions  are  allowed  to  vibrate,  so  the  gate  is  essentially 
local  and  the  gate  complexity  does  not  depend  on  the  crystal 
size,  (c)  The  relative  response  of  the  ions  for  the  gate  shown 
in  (b)  (characterized  by  the  largest  spin- dependent  shift  \q^\ 
during  the  gate  time  r).  We  take  a  relative  unit  where  \q^\ 
for  the  target  ions  haven  been  normalized  to  1.  The  fast  decay 
of  the  response  as  one  moves  away  from  the  target  ions  (59th 
and  62nd)  shows  that  the  gate  involves  vibration  of  only  local 
ions. 


U{t) 


exp 


+  an*('rHK 


+  i  F  ())mn{T)<JZm<JZn  , 

m<n 


(2) 


where  a^(r)  =  /q  Xn(^)ffnelUJkt dt  characterizes  the  resid¬ 
ual  entanglement  between  ion  n  and  phonon  mode  k  and 
<?W(t)  =  2/0Tdt2/o2  dti  Efefi'm5nXm(i2)Xn(ii)sinwfc(i2  - 
ti)  represents  the  effective  qubit-qubit  interaction  between 
ions  m  and  n. 

For  a  two-qubit  gate  on  an  ion  pair  i  and  j,  we 
direct  laser  light  exclusively  on  these  two  ions  (f^(t)  = 
Qj(t)  =  Q(t)  and  all  other  f}n(t)  =  0),  and  the  evolution 
operator  reduces  to  the  standard  controlled  7r-phase  (CP) 
gate  for  a^(r)  =  0  and  (j)jn(r)  =  7t/4.  For  a  large  ion 
crystal,  the  residual  entanglement  with  the  motional 
modes  cannot  be  eliminated  completely,  but  we  can 


minimize  the  resulting  gate  infidelity  by  optimizing  the 
laser  pulse  shape  Q(t)  [23].  Assuming  each  phonon  mode 
k  is  cooled  to  temperature  Z^,  the  infidelity  of  the 
CP  gate  from  the  residual  motional  entanglement  is 
given  by  5F  =  [6  —  2(r$  +  Tj)  —  T+  —  F_]/8  [23],  where 

U0)=exp[-Efe  lafo-)(r)l2/W2L  r±=exP[-Efc  \ai(T)± 

aj(r) |2/3fc/2],  and  (3k  =  cofh{hu>k/kBTk). 

To  minimize  the  gate  infidelity  SF ,  we  break  the  laser 
pulse  on  the  two  ions  into  uniform  segments  of  constant 
intensity  as  shown  in  fig.  3(b)  and  optimize  the  values 
(i  =  1, . . . ,  M)  over  M  equal-time  segments  [23].  The 
control  of  Q(t)  from  one  segment  to  the  next  can  easily  be 
accomplished  with  optical  modulators.  (Alternatively,  we 
can  modulate  the  detuning  //  of  the  laser  pulse  as  a  control 
parameter.)  After  optimization  of  fiW,  the  infidelity  SF 
is  shown  as  a  function  of  the  detuning  /jl  in  fig.  3(a)  for 
M  =  5  segments.  For  this  example,  we  perform  a  CP  gate 
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on  the  59th  and  62nd  ions  in  this  120-ion  chain.  With 
an  appropriate  choice  of  fiW  and  /i,  the  infidelity  can  be 
made  negligible  (well  below  10-5).  For  this  calculation,  we 
assume  Doppler  cooling  for  all  modes  and  take  gate  times 
r  in  the  range  50ro  to  500ro,  where  To  =  27t/ujx  =  0.2 /xs  is 
the  period  of  transverse  harmonic  motion.  The  gate  can 
certainly  be  faster  with  stronger  laser  beams  (there  is  no 
speed  limit),  and  with  a  faster  gate,  the  control  becomes 
easier  as  the  gate  becomes  more  localized  (fig.  3). 

Interestingly,  we  use  only  a  few  control  parameters 
(M  =  5  segments)  to  perform  a  high-fidelity  gate  that 
involves  excitation  of  hundreds  of  transverse  phonon 
normal  modes.  This  is  possible  because  the  gate  has  a 
local  character  where  the  contribution  to  the  CP  gate 
comes  primarily  from  the  spin-dependent  oscillations  of 
the  ions  that  are  close  to  the  target  ions.  To  show  this,  we 
plot  the  response  of  each  ion  in  fig.  3(c)  during  the  gate 
operation.  Note  that  the  displacement  qn  of  the  n-th  ion  is 
spin-dependent  during  the  gate,  and  we  can  use  its  largest 
magnitude  \q^\  over  the  gate  time  r  to  characterize  the 
response  of  ion  n,  as  is  shown  in  fig.  3(c).  The  ion  response 
decays  very  fast  from  the  target  ions  (59th  and  62nd  in 
this  case)  and  can  be  safely  neglected  after  a  distance  of 
a  few  ions.  Thus  during  a  gate,  only  the  motion  of  ions 
near  the  target  ions  are  important,  and  the  other  ions 
largely  remain  in  their  equilibrium  positions.  The  resultant 
control  parameters  from  this  approximation  are  almost 
identical  to  those  shown  in  fig.  3(b).  Owing  to  the  local 
character  of  the  gate,  the  complexity  of  a  gate  operation 
does  not  depend  on  the  chain  size,  and  we  can  perform 
gates  in  parallel  on  ions  in  different  regions  of  a  large  chain. 

As  we  use  primarily  the  quasi-local  phonon  modes  for 
the  gate  operations,  we  assume  only  local  entangling 
gates  where  the  distance  between  the  target  ions  is 
small  compared  with  the  length  of  the  whole  ion  chain. 
(For  distant  quantum  computing,  gate  operations  can 
be  implemented  through  aids  of  a  series  of  mediate 
local  SWAP  gates.)  Recent  studies  have  shown  that 
with  only  nearest-neighbor  entangling  gates  in  a  two- 
dimensional  (2D)  lattice,  the  error  threshold  for  fault- 
tolerant  computation  can  be  still  very  good,  close  to  the 
level  of  one  percent  [24].  In  a  one  dimensional  (ID)  ion 
chain,  if  we  perform  gates  with  distance  up  to  (which 
is  still  relatively  local  compared  with  the  size  N  of  the  ion 
chain),  this  simulates  a  2D  system,  and  the  error  threshold 
for  this  case  should  be  as  least  as  good  as  the  2D  case  with 
nearest  neighbor  entangling  gates.  Moreover,  in  the  ID 
case,  it  is  even  possible  to  perform  fault-tolerant  quantum 
computing  with  only  next  to  nearest-neighbor  entangling 
gates  [25],  albeit  with  a  more  demanding  threshold. 

We  now  discuss  several  sources  of  noise  for  gates 
in  a  large  ion  crystal  and  show  that  their  effects  are 
negligible.  First,  the  axial  ion  modes  have  large  phonon 
occupation  numbers  under  Doppler  cooling  alone,  and 
the  resulting  thermal  spread  in  position  along  the  axial 
direction  can  degrade  the  effective  laser  interaction.  For 
example,  the  lowest  axial  mode  in  a  120-ion  chain  of  Yb+ 


ions  with  a  spacing  z  ~  10  fim  has  a  frequency  of  only 
<^lo/2tt  =  9.8  kHz  and  a  mean  thermal  phonon  number 
no  ~  7/^lo  ~  103  under  Doppler  laser  cooling  (radiative 
half-linewidth  7/27T  =  10 MHz).  We  assume  the  quantum 
gate  laser  beams  are  directed  along  the  transverse  direc¬ 
tion  with  an  axial  Gaussian  laser  profile  f 1(z)  oc 
centered  on  each  ion.  The  beam  waist  is  taken  as 
w  =  z/2.5  «  4 /mi  so  that  the  cross-talk  error  probability 
between  adjacent  ions  is  Pc  =  <  10-5.  The  posi¬ 

tion  fluctuation  5zn  of  the  n-th  ion  causes  the  effective 
Rabi  frequency  to  fluctuate,  resulting  in  a  gate  infidelity 
SFi  «  (tt2 /4:)(5Qn/Qn)2  «  (it2 /4)(8zn/w)A .  The  fluctua¬ 
tion  8zn  can  be  calculated  exactly  from  summation  of 
contributions  of  all  the  axial  modes  and  its  value  is  almost 
independent  of  the  index  n  for  the  computational  ions  (see 
Appendix).  Under  Doppler  laser-cooling,  Szn~  0.26 //m 
and  the  corresponding  infidelity  is  SFi  =4.4  x  10-5.  The 
position  fluctuation  of  the  ions  may  also  lead  to  anhar- 
monic  ion  motion,  whose  contribution  to  the  gate  infidelity 
can  be  estimated  by  8F2  (W*)2~  6-8  x  10  4.  Finally, 

in  the  transverse  direction  we  estimate  the  infidelity 
caused  by  the  higher-order  expansions  in  the  Lamb-Dicke 
parameter.  As  all  the  transverse  modes  have  roughly  the 
same  frequency  uik  ~  the  effective  Lamb-Dicke  para¬ 
meter  for  the  transverse  modes  is  rjx  =  |  Ak|  ^h/2mu)x  « 
0.038  for  Yb+  ions  at  ux/2tt  =  5  MHz,  with  each  mode 
containing  a  mean  thermal  phonon  number  nx~  2.0  under 
Doppler  cooling.  The  resultant  gate  infidelity  is  estimated 
to  be  SFs^  7r27y4  (n2  +  nx  +  1/8)  «  7  x  10-4  [19,26] .  Note 
finally  that  sideband  cooling  is  possible  in  the  transverse 
direction  as  all  the  modes  have  nearly  the  same  frequency, 
thus  reducing  the  gate  infidelity  due  to  transverse  thermal 
motion  by  another  order  of  magnitude. 

In  summary,  we  have  shown  through  explicit  examples 
and  calculations  that  it  is  feasible  to  stabilize  large 
linear  ion  crystals  where  the  gate  complexity  does  not 
increase  with  the  size  of  the  crystal  and  the  gate  infidelity 
from  thermal  fluctuations  can  be  made  negligibly  small 
under  routine  Doppler  cooling.  The  results  suggest  a 
realistic  prospect  for  realization  of  large-scale  quantum 
computation  in  a  simple  linear  ion  architecture. 
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Appendix:  thermal  position  fluctuation  of  the 
ions  along  the  axial  direction.  —  We  address  indi¬ 
vidual  ions  through  focused  laser  beams  which  typically 
take  a  Gaussian  shape  along  the  ^-direction  with  £Ln{z)  oc 
e~zn  lw  5  where  zn  =  z  —  zn  is  centered  at  the  equilib¬ 
rium  position  zn  of  the  n-th  ion.  Under  the  Doppler 
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Fig.  4:  The  axial  position  fluctuation  5zn  is  plotted  along 
the  ion  chain,  which  is  about  0.26 /xm  (averaged)  for  the 
computational  ions  (n  from  11  to  110). 


temperature,  the  ions  have  significant  thermal  fluctua¬ 
tion  of  their  positions  along  the  ^-direction,  which  leads 
to  an  effectively  fluctuating  laser  amplitude  Qn(z)  and 
induces  infidelity  to  the  gate  operation.  This  position 
fluctuation  influences  both  the  single-bit  and  the  two- 
bit  operations  in  the  same  way.  To  quantify  the  gate 
error  caused  by  this  fluctuation,  let  us  consider  a  spin- 
flip  gate  operated  on  the  n-th  ion  as  a  typical  exam¬ 
ple.  For  a  spin-flip  with  a  7r-pulse,  the  gate  fidelity  is 
given  by  Fi  =  sin2(finr  +  5Q.ht)  «  1  -  (7r2/4)(<5f2n/Cra)2, 
where  Qn  is  the  expectation  value  of  the  Rabi  frequency 
which  satisfies  Dnr  =  7r/2  for  a  spin-flip  gate,  and 
is  its  fluctuation  caused  by  the  position  fluctuation  of 
the  ion.  From  Cln(z)  oce~z™/w  «1  —  z'2/w2  around  the 
equilibrium  position,  the  gate  infidelity  SF±  =  1  —  F\  = 
(n2  /  A)  (6  zn  /  w)4 ,  where  5zn  =  (z'4  —  )1^4  characterizes 

the  thermal  position  fluctuation  of  the  n-th  ion  along  the 
axial  direction. 

The  phonon  modes  are  in  thermal  equilibrium 
under  the  Doppler  temperature  T,  with  their  density 
operator  given  by  pm  =  ]\kY;{nk}  Pk\nk)(nk\,  where 
Pk  =  nkk  /(rife  +  l)n/c+1  is  the  probability  of  having 
nk  phonons  in  the  fc-th  mode,  and  nk  =  / (hwk) 

is  the  average  phonon  number.  From  z!£  =  tr(z'2pm) 
and  z'4  =  tr(z^pm),  we  explicitly  have  5zn  = 

^/h/2m[\/2  Yjk  (K’k)2  (2nfe  +  l)/w2ifc]1/2  where  ojz,k 
and  b^k  denote  the  eigenfrequencies  and  eigenmatrices 
of  the  axial  modes.  For  our  example  of  a  120- ion  chain 
with  the  Doppler  temperature  &#T/ft  =  62MHz  for 
the  Yb+  ions,  the  resultant  position  fluctuation  5zn  is 
plotted  in  fig.  4  for  all  the  ions.  One  can  see  that  for  the 


computational  ions  (n  from  11  to  110),  Szn~  0.26 //m 
with  its  value  almost  independent  of  the  ion  index.  The 
position  fluctuation  is  still  significantly  smaller  than  the 
ion  spacing  («  10  gm),  which  ensures  a  tiny  gate  infidelity 
SFi  as  discussed  in  the  main  text. 
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